UNCLASSIFIED 


.(  405  826 


DEFENSE  DOCUMENTATION  CENTER 

FOR 

SCIENTIFIC  AND  TECHNICAL  INFORMATION 

CAMERON  STATION.  ALEXANDRIA,  VIRGINIA 


UNCLASSIFIED 


NOTICE:  When  govenment  or  other  drawings,  speci¬ 
fications  or  other  data  are  used  for  any  puipose 
other  than  In  connection  with  a  definitely  related 
govemnent  procurement  operation,  the  U.  S. 
Government  thereby  incurs  no  responsibility,  nor  any 
obligation  irtiat soever;  and  the  fact  that  the  Govern¬ 
ment  may  have  fozmilated,  furnished,  or  In  any  way 
supplied  the  said  drawings,  specifications,  or  other 
data  Is  not  to  be  regarded  by  Implication  or  other¬ 
wise  as  In  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rl^ts 
or  permission  to  manufacture,  use  or  sell  any 
patented  Invention  that  may  In  any  way  be  related 
thereto. 


RADC‘TDR-63-124 


MAY  1963 


00 

o 


DISCRETE,  FINITE  REPRESENTATION  OF 
A  LINEAR,  STATIONARY  SYSTEM 
By 

P.  J.  Butterly 


THE  JOHNS  HOPKINS  UNIVERSITY 
CARLYLE  BARTON  LABORATORY 
Charles  and  34th  Streets 
Baltimore  18,  Moryland 


TECHNICAL  DOCUMENTARY  REPORT  NO.  RADC-TDR-63-124 


Contract  No.  AF  30(602)-2597 
Program  Element  Code -62405454 -7600 


1 1  ■  I '  i  i~>J  i  I  i\  I.  ■ 


Prepared  For 

Rome  Air  Development  Center 
Air  Force  Syitemi  Conmond 
United  States  Air  Force 
Griffisi  Air  Force  lose 
New  York 


PATENT  NOTICE;  When  Government  drawings,  specifications, 
or  other  data  are  used  for  any  purpose  other  than  in  connection  with 
a  definitely  related  Government  procurement  operation,  the  United 
States  Government  thereby  incurs  no  responsibility  nor  any  obliga¬ 
tion  whatsoever  and  the  fact  that  the  Government  may  have  formulated, 
furnished,  or  in  any  way  supplied  the  said  drawings,  specifications 
or  other  data  is  not  to  be  regarded  by  implication  or  otherwise  as  in 
any  manner  licensing  the  holder  or  any  other  person  or  corporation, 
or  conveying  any  rights  or  permission  to  manufacture,  use,  or  sell 
any  patented  invention  that  may  in  any  way  be  related  thereto. 


Qualified  requesters  may  obtain  copies  from  DDC,  Arlington  Hall 
Station,  Arlington  12,  Virginia.  Orders  will  be  expedited  if 
placed  through  the  librarian  or  other  person  designated  to  request 
documents  from  DDC. 


RADC-TDR-63-124 


MAY  1963 


DISCRETE,  FINITE  REPRESENTATION  OF  A 
LINEAR,  STATIONARY  SYSTEM 

by 

P.  J.  Butterly 


THE  JOHNS  HOPKINS  UNIVERSITY 
CARLYLE  BARTON  LABORATORY 
Charles  and  34th  Streets 
Baltimore  18,  Maryland 


TECHNICAL  DOCUMENTARY  REPORT  NO.  RADC-TDR-63-124 


Contract  No.  AF  30(602)-2597 
Program  Element  Code:  62405454 
Project  Number  4505 
Task  Number  450501 


Prepared  for 


Rome  Air  Development  Center 
Air  Force  Systems  Command 
United  States  Air  Force 

Griffis s  Air  Force  Base 
New  York 


PATENT  NOTICE;  When  Government  drawings,  specifications, 
or  other  data  are  used  for  any  purpose  other  than  in  connection  with 
a  definitely  related  Government  procurement  operation,  the  United 
States  Government  thereby  incurs  no  responsibility  nor  any  obliga¬ 
tion  whatsoever, and  the  fact  that  the  Government  may  have  formulated, 
furnished,  or  in  any  way  supplied  the  said  drawings,  specifications 
or  other  data  is  not  to  be  regarded  by  implication  or  otherwise  as  in 
any  manner  licensing  the  holder  or  any  other  person  or  corporation, 
or  conveying  any  rights  or  permission  to  manufacture,  use,  or  sell 
any  patented  invention  that  may  in  any  way  be  related  thereto. 


FOREWORD 


This  report  was  prepared  by  Carlyle  Barton  Laboratory,  The 
Johns  Hopkins  University,  Baltimore,  Maryland  on  Air  Force 
Contract  No.  AF  30(602)_2597  under  Project  No.  4505  of  Task  No. 
450501.  The  work  was  administered  by  the  Electronic  Warfare 
Laboratory,  Rome  Air  Development  Center.  Mr.  Haywood  E, 

Webb  was  the  project  engineer. 


PUBLICATION  REVIEW 


CO 


This  report  has  been  reviewed  and  is  approved. 


Approved; 


EDWARD  N.  MUNZER  } 

Chief,  Electronic  Waiuare  laboratory 
Directorate  of  Intelligence  &  Electronic  Warfare 


Approved: 


ROBERT  M  QUINN,  J^.,  Col,  USAF 
Directoijf'^of  Intelligence  &  Electronic  Warfare 


-iii- 


‘1 


ABSTRACT 

This  report  describes  a  linear  system  representation  suitable  for 
use  with  signals  of  finite  dimensionality.  This  representation  takes 
the  form  of  a  matrix  with  elements  which  depend  both  on  the  choice  of 
signal  basis  and  on  the  transmission  properties  of  the  system.  The 
report  is  divided  into  two  parts.  The  first  part  is  devoted  to  the  deri- 
CNJ  vation  of  the  system  matrix  for  an  arbitrary  system  transfer  function 

H 

for  several  widely-used  signal  bases.  The  second  part  consider  s  the 
^  case  in  which  only  the  system  matrix  is  known.  An  equivalent  repre¬ 

sentation,  expressedin terms  of  a  continuous  parameter,  is  introduced 
to  facilitate  approximation  of  the  system  transfer  function,  and  trans  - 
formation  formulas  are  established  for  a  number  of  cases  where  they 
have  a  simple  form. 
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1.  DISCRETE,  FINITE  REPRESENTATION  OF  A  LINEAR 
STATIONARY  SYSTEM 
1.1.  INTRODUCTION 

The  efficient  representation  of  a  given  signal  class  in¬ 
volves  the  selection  of  a  minimal  set  of  component  functions  in  terms 
of  which  every  member  of  the  class  may  be  described  with  acceptable 
accuracy.  In  many  instances,  as,  for  example,  in  the  representation 
of  .signals  limited  in  time  and  frequency,  a  finite  set  of  components  is 
sufficient  for  the  complete  representation  of  the  class.  In  order  to 
devise  a  characterization  of  a  linear,  stationary  transmission  system 
compatible  with  this  method  of  signal  representation,  consider  the  re¬ 
sponse  of  this  system  to  an  arbitrary  signal  chosen  from  such  a  signal 
class.  In  general,  it  will  not  be  possible  to  obtain  a  complete  represen¬ 
tation  of  the  response  in  terms  of  the  original  set  of  component  functions; 
however,  the  least-mean-square  approximation  of  the  response  which 
may  be  so  represented  leads  to  a  satisfactory  characterization  of  the 
system  provided  the  resulting  error  is  small  enough  .  Where  it  is  not, 
a  re-examination  of  the  set  of  signal  components  chosen  becomes 
necessary  and  a  larger  set  may  have  to  be  used. 

Clearly  this  situation  corresponds  to  the  notion  of  a  linear 
transformation  of  the  original  signal  space  on  to  itself  and  to  the  re¬ 
presentation  of  the  transmission  properties  of  the  system  by  the  matrix 

2 

of  the  transformation  .  The  elements  of  this  matrix  will  depend  both 

on  the  choice  of  basis  -  that  is,  on  the  signal  components,  and  on  the 

system  itself.  The  derivation  and  manipulation  of  this  matrix  is  the 

subject  of  Part  1  of  this  report.  The  treatment  parallels  that  given 

3 

in  a  recent  Internal  Memorandum  . 

Let  |f)  be  an  arbitrary  signal  chosen  from  the  given 
signal  class.  Then  |F)  maybe  completely  represented  in  terms  of 
some  basis  set 

k=l„2, . N 


I'  ; 


I 


I 
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N 

that  is,  |F)  =  ^ 

k=l 

If  |h|  is  the  linear  operator  referred  to  above, then,  following 
1 


Huggins  ,  we  can  write 
N 

I  v'v 

j  =  l 


(1.1) 


The  left  side  of  this  relationship  is  the  true  response  of  the  system. 

The  right  side  is  a  vector  lying  in  the  original  signal  space.  Clearly 
2 

the  N  values  of  =  1.  2,  ....  N,  for  which  the  equivalence 

is  in  the  least-squares  sense,  constitute  a  discrete  representation 
of  |h|  which  provides  a  least-squares  approximation  of  the  system 
response.  In  order  to  ascertain  those  values  of  h^j^  let 

N 

l'k>  =  -  I 

j=  1 

then,  for  minimum  mean-square  error:- 

j=i  j=i 

N  . 


^  <~fi®k>=  -  <^pi[iHi\>  -  y 

pq  j=l 


N 


=  0  when  <^p|H!\>  =  Y 

j=  1 


m 


124 


1^4 
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m  V  =  iHk.  > 

g  .  =  (ip  |4.> 
t-pj  ^  '  4 

For  minimum  error, 

=  X  V«Pi 

j=  1 


m 


pk 


or,  in  matrix  notation: - 
M  =  G  H 


and 


H  =  G  ^  M 


(1.2) 


(1.3) 


H  is  the  required  system  matrix  with  elements  h^j^  j  =  1,  2,  .  .  .N  and 
G  is  the  familiar  Gram  Matrix. 

As  an  alternative  procedure,  an  orthonormal  basis  could  have 
been  derived  by  means  of  a  linear  transformation  of  the  given  basis, 
in  which  case  the  Gram  Matrix  would  reduce  to  the  Unit  Matrix.  This 
suggests  a  relation  between  the  Gram  Matrix  and  the  matrix  which 
transforms  a  given  basis  to  an  orthonormal  one.  This  relationship 
is  explored  in  Appendix  1. 


To  determine  the  Error  Energy: - 
N 


<~(i'k>=  -  I  4  (^ii]  [i«i+k>-  I 


i=  1 
N  N 


N 

j=  1 


m5<|h|h|>|.^)  +  ll  huhjkSij 

i=  1  j=  1 


I 


-1  V  - 1  ''h  ”ik 

j=  1  i  =  1 

N  N 


i=  li=  1 


-I  "u 


ii  “ik 


where  |r,  )=  |h|  I,  > 


N  N 


=  -  I  I 


i=li=l 


>!<  -1 

1..  g  ..  m., 

ji  ®  ji  ik 


E=U-M^  where  ^  I  ^k^ 


N  N 


Then  the  error  energy  "  ^  ^  ^  ®k^  j'^st  the  Bum  of  the 


elements  of  E, 


f  =  1  k=l 


For  an  orthonormal  basis,  G  =  1,  M  =  H  in  which  case:- 
N  Tl 

Error  Energy  =  ^  ^~k^^k^  ”  Z  ^  h^, which  is  the  true 

k=l  k=lj=l 


output  energy  less  the  sum  of  the  squares  of  the  matrix  elements. 


I 
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1.2  LINEAR  TRANSFORMATIONS  AND  THE  DIAGONAL 
REPRESENTATION 

It  is  often  possible  to  simplify  the  system  matrix  by  means 
of  a  linear  transformation  of  the  signal  basis. 

Consider  any  two  complete  bases  and  i  =  1,  2,  .  ,  .N 

and  let 


but 


and 


If 


IV  = 

N 

I  'jkl+j' 

|F>  = 

I'.. 

N 

k=  1 

N 

=  I 

4  I 

k=l 

i’i 

N  N 

=  I  I  'jk4i+j>' 

j  =  i  k=i 

|F)=  ^  , 

j=  1 
N 

4^  =  Y'  „  A 


J  L  jk  k 
k=l 


R>4,  =  H4,£)  + 


then 


1 


But 

5.>^,  =  Di 

t4  =  c-‘H^£  ■  (2.1) 

If  can  be  diagonalized,  say  by  the  similarity  transformation 

=  ,  (2.2) 

then,  from  Equations  (2. 1)  and  (2.  2),  it  is  clear  that  C  =  ^  gives  the 
linear  transformation  of  the  signal  basis  which  will  always  lead  to  a 
diagonal  form  of  the  system  matrix. 


I'i 

CO,! 


! 

( 

i: 

i 


V  ,! 
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1.3  EXPONENTIAL  BASES 

The  choice  of  basis  is  determined  primarily  by  the  nature  of 
the  input  signal  which  must  be  adequately  represented.  So  far, 
error  estimation  has  been  on  the  assumption  that  the  description  of 
the  input  vector  is  complete.  For  an  extremely  wide  class  of 
signals,  however,  growing  (or  decaying)  exponentials  provide  a 
suitable  characterization.  Apart  from  the  discontinuity  resulting 
from  the  signal  epoch,  these  components  are  eigenfunctions  of  the 
system.  As  will  be  seen,  this  leads  to  a  diagonal  system  matrix 
with  easily  computed  elements  provided  the  exponents  are  distinct. 
The  case  of  a  characterization  based  on  a  single  exponent  is  dis¬ 
cussed  separately.  A  more  general  discussion  of  the  use  of 

1  4 

exponential  components  will  be  found  in  Huggins  and  Lai  . 

A.  Forward  Components 
s- t 

Let  v|jj^(t)  =  e  ^  t  >  0,  Sj^  =  Cj^  +  j(Dj^ 

4'k(s)  =  0  £  a  s  =  a  +  juu 

k 

”^fk“  J  ^  s-s  Zirj  • 

p,  ~  S  *"  S  ^  K 


From  Appendix  Z,  F  may  be  taken  as  the  juj  axis 
H(-s^) 

m.,  =  — jj; -  for  H(s)  having  no  poles  in  right  half  plane. 

+  ^k 


Also, 


Sfk 


"f+®k 


In  matrix  notation: - 


M  =  Hg  G 


where 
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H  = 
— s 


H(-s 


0 


0 


0  0  - 

H(-S2)  0  - 

0  H(-s*)  - 


and 


From  Equation  (1.2) 

H,  =  G‘^  H  G 
—  — s  — 


H  =  G  H,  G' 
— s ^  — 


(3.1) 

(3.2) 


From  (2.  2),  the  change  of  basis  which  will  result  in  a  diagonal 
system  matrix  is  given  by;- 


B.  Backward  Components 
Let  =  e 

*k<’»  =  • 

k 


t  <  0  >  0 

as  0 


=  H(s)  ^ 

r  s+s , 


ds 

-s+Sr  ^Trj 


H(s, 


Sk+S| 


r  s+s 


1  ds 
^  -s+s^_  2lTj 


»!< 
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M  =  G  H 
- - 8 

where 


H  = 
— s 


H(s^) 

0 

0 

0 


0 

H(s^) 

0 

0 


0 

0 

H(83) 

0 


o:( 


and 

H,  =  G"^  M  =  H  ,  (3,  3) 

showing,  that  for  backward  components  this  choice  of  basis  leads 
directly  to  a  diagonal  system  matrix. 


CO 


V 


C.  Matrix  Elements  for  a  Basis  Consisting  of  Laguerre  Functions 


The  reduction  of  the  system  matrix  to  a  diagonal  form 

results  from  the  choice  of  distinct  exponents  in  the  basis.  There  is 

one  important  case  in  which  this  does  not  occur,  namely,  when  the 

input  signal  is  characterized  in  terms  of  a  single  exponent.  In  this 

1  5 

instance,  the  orthonormal  basis  derived  by  the  Kautz  method  ’ 
has  a  particularly  simple  form.  Consider  the  case  of  Sj^  =  a,k=l,  2, 
.  .  .  .  N  (forward  components) 


=  (-a  -q 


1/2 


(s+a  ) 


*,k-l 


(s  -a 


k  =  1,2 . N 


"fk 


(-s+a) 


f-1 


I  j.  ’\k-l 
(s+a  ) 


=  (-a -a  )  r  H(s)  - - ^ 


r  (-8-a  ) 


(s-a)' 


ds 

2itJ 


-10 


=  -(-a-a  )  f 


H(s)(s-a)^  ^  ds 


r  (s+a  ) 


FTI  2iTj 


where  6  =  i  -  k 


-( -a  -a  ) 

8s 


s  =  -a 


=  0 


6  >  0 
6  <  0 


(3.4) 


h  =  H(-a  ) 
o  '  ' 

h.  =  (-a-a*)  HM-a"^) 


i,  -  -Cl 

^2 - Z— 


2H'(-a*)  +  (-a-a*)  H"(-a*)] 


_  (-g-g  ) 


a^r,  ’i‘2  *“1 

^'[6H'{-a  )  +  6(a-a  )H"(-a  )  +  (-a-a  )  H"'(-a  )J 


etc. 


The  H  matrix  is  therefore  triangular  and  the  values  of  its  elements 
depend  only  on  their  distance  from  the  diagonal 

If  a  =  -1/2 

(3.5) 

s  =  -l/Z 

6  i  0 

=  0  6  <  0 


[(.MA)‘-‘  H(.)] 
ds 


h^  =  H(l/2) 
h^  =  H'(l/2) 

h^  =(l/2)  2H'(l/2)+H"(l/2)] 

h3=(l/6)  6H''(1/2)+6H"(1/2)+H"'(1/2)]  etc. 


03 
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1.4,  MATRIX  ELEMENTS  FOR  A  BASIS  CONSISTING  OF 
ELEMENTARY  "HOLD"  FUNCTIONS 

Exponential  functions  are  suitable  for  characterizing  signals  of 
semi-infinite  duration  especially  when  there  is  prior  knowledge  of  the 
distribution  of  complex  frequencies.  They  do  not  lend  themselves  to 
an  exact  description  of  a  time -limited  signal  such  as  might  result 
from  the  application  of  quantization  or  pulsed- code  techniques.  The 
introduction  of  compound  errors  may  be  avoided  by  using  the  so-called 
"B"  basis,  defined  as  follows:- 

’  f 


4Jj^(t)  =  1  (k-l)st<k 

=  0  t  <  (k-1),  t  2:  k 

k  =  1,2, _ N.. 


This  is  an  orthonormal  basis.  Considering  forward  components  only:- 

,1,  ,  ,  _  1  ,,  -s.  -{k-l)s  ^ 

+j^(s)  -  -  (l-e  )  e  '  -oo  a 


II 

.i: 

a 

|h1\> 

J  l-e« 

e^j-^)®H{s) 

r 

r  H(s)j 

f  s  .  -s  , 

l^e  +  e  -  2 

-e  ®  -(k-l)s  ds 

~s  ®  ZttJ 


6b  ds 
®  2irj  ’ 


(4.1) 


where  6  =  j-k,  and  F  may  be  taken  as  slightly  to  the  right  of  the  ju) 
axis.  Again,  the  values  of  the  elements  depend  only  on  distance 
from  the  diagonal. 

Case  1  6  <  0 

hg  =  0  since  the  integrand  of  (4. 1)  is  analytic  to  the 

right  of  r. 


I 


Case  Z 


6  =  0,  diagonal  terms 


=  I 


H(s)  s  ds 

—7~  ®  2irj 
r  s  •' 


he  =  ^^(0)  +  ^^(1)  ^\o),  by  Appendix  3. 


(4.2) 


Case  3 


6  >  0 


r  s 


P  H(s)  (6+l)s  ds  ,  p  H(s)  6b  ds  ,  p  H(s)  (6-l)s  ds 

~  ~T-  ^  2^J"^J  “T-®  ZSJ 

rs  *'rs  “'rs  •’ 


Using  the  result  in  Appendix  3,  we  have  that 

hg  =  j’-(6-l)h^'^^(0)  +  h^’^\6-l)-h^‘^\0)]- 
-2  |-6h^‘^)(0)  +  +  h^'^^(O)]- 

+  ■[-(6+4)h^'^^(0)  +  h^'^^(6+l)-h^"^\o)]- 

h,  =  h^'^\6+l)  -  2h^‘^‘\6)  +  h^‘^^(6-l) 


(4.  3) 


This  completes  the  determination  of  matrix  elements  for  signals 
existing  only  for  positive  time.  Otherwise,  this  matrix  is  one  of  the 
four  sub -matrices  necessary  to  describe  the  system.  For  the  other 
three. 


hj^=  (^Pj^IhK.)  =hj^. 

i—®  '  S  ZlTJ 

.  j  hjj^  =  0  for  all  j  and  k 


(4.4) 


(4.5) 


1 


-13. 


"jk  -s 


l-e 


1-e®  (k-l)s  _d£_ 
-  s  ®  Zirj 


r  s 


J 

H(s)  fg(6-2)s 

I  i  2irj 


6  =  j  -  k  =  j  +  k 


6  s  2 


hj^={-(6-2)h^"^\o)  +  h^'^^(6-2)  -  h^'^^IO)} 

-2  j-(6-l)h^‘^^(0)  +  h^'^^(6-l)  -  h^'^^(O)]- 

f 

it 

ii 

j-6h^'^\o)  +  h^"^^(6)  -  h^"^^(0)]- 

i', 

i: 

CM  ^  ^ 

=  h^~^^{6)  -  h^’^^(6-l)  +  h^'^^(6-2) 

(4,6) 

fe  ■ 

CO 

6  :a  2 

A.  Effect  of  Mid-Interval  Sampling 

,  f 

) 

j 

It  is  often  more  convenient  to  use  the  basis  defined  by:- 
i|jj^(t)  =1  (k  -  1/2)  i  t  <  (k  +  1/2) 

=  0  t  <  (k  -  1/2),  t  ^  (k  +  1/2) 

k  =  0,  1,  2,  .  .  .  .N 


A  8 

l-e 


'  8 

6s  r  8  -  ,  -s'!  ds 

e  -2+e  I  JJ-. 


ds 

2Trj 

(4.7) 


6  =  j  -  k 
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Thus,  for  a  given  6,  these  matrix  elements  are  the  same  as  those 
determined  in  Section  1. 4,  However,  both  j  and  k  can  now  be  zero. 
As  before, 

^jk  "  \j  6=7-k  =  k-  j) 

hjj^  =  0  all  j,  k,  except  j  =  k  =  0.  6  =  j  -  k  =  -(j  +  k)  <  0. 

,  -  f  1  /i  Sv  ~^\  (k+l/2)B  ds 

h.c- -  — (1-e  )e  H(s)^l-e  )e  '  -tt— r 

jk  -s  '  '  '  s'  '  2irj 

=  L  -T-  -2+e  }e  ^  (4.8) 

X  s 

where 

6  =  j  -  k 
=  j  +  k 

h.r  =  h^‘^^(6+l)  -  2h^'^\6)  +  h^'^^(6-l)  for  6  i  2 

from  Equation (4.  3). 


Summarizing,  we  have 


hg  =  0,  6  <  0 

hg  =  -h^'^^(O)  +  h^'^^(l)  -  h^"^^{0)  ,  6=0 

hg  =  h^"^\6  +  1)  -  2h^"^^6)+h^'^^(6-l),  6  >  0  , 

where  6  is  defined  as  above  for  each  submatrix. 

Only  in  elements  such  as  h^^  does  the  system  matrix  differ 
from  that  of  Section  1, 4,  although  it  differs  in  size  by  at  least  one 
row  and  column.  In  integrating  the  impulse  response  for  insertion 
in  the  above  formulas  arbitrary  constants  are  introduced  but 
do  not  appear  in  the  final  evaluation  of  the  matrix  elements. 
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1.5  MATRIX  ELEMENTS  FOR  A  BASIS  CONSISTING  OF 
ELEMENTARY  TRIANGULAR  FUNCTIONS 


A  complete  description  of  any  function  which  can  be  reduced 
to  straight  line  segments  is  possible  using  the  elementary  triangular 
component  defined  by:- 

^^{b)  =  -  Z  +  e~^)  e~^^,  k  =  0,  1,2,  ...N. 

s 

A.  The  M  Matrix 


-2+e®)ej®H(8) 


(e®-2+e 

s 


de 

lirj 


=  J  H{s2g8{6+2)  _  ^^s(6+l)|^,,6s  ^^8(6-1)^  ^8(6-Z)~|dB 
r  s”*  L 


(5.1) 


Case  1 


6  <  -1 


plane 
Case  2 


m., 


0  since  the  integrand  is  analytic  in  the  right  half 


6  =  -1 


m 


jk 


H(s) 

~T~ 


3 


e 


ds 

Ztrj 


1 

F 


h^'^\o)  - 


by  Appendix  3 . 
Case  3  6=0 


From  Equation  (5.1) 


m 


6 


ds 

2irj 


=  -  |.h^'^)(:0)  -i.h^‘^)(0)  -  2h^'^^{0)  -  h^’^)(0)+h^'‘^\2) 


I 


.4j^-  -  h^'^^{0)  -  h^"'^^(0)  +  h^''^^(l)] 

=  _|.h^"^)(0)+  2h^"^\o)+  3h^"^\o)  +  h^'^^Z)  -  4h^''^^(l)  , 


Case  4 


6=1 


M<|)|e^".4e^%  6e®}  ^ 


=  h^'^\o)  -  |-h^"^^(0)  -  3h^‘^\o)-h^"'*^(0)+h^‘'^^(3) 

-4[-  |h^  ~^\q)  - 1  h^  '^^0)  -2h^  V)-h^  '^^(0)+h^  ■‘^^(2)] 

+6[-  -l-h^  ‘^^(0)  -  ih^  ■^^(O)-h^  ■^^(O)-h^  "^^(0)+h^ 


=  .^h^"^){0)  +  |h^"^)(0)-h^'^^(0)-3h^"^V)+h^'‘^^(3) 


-4h^‘^^(2)  +  6h^'^^(l) 


Case  5 


From  Equation  (5,1) 

mg=-i5-^^h^'^\o)-^-^^^h^"^^(0)-{6+2)h^'^^(0)-h^"'^^(0)+h^'^^(6+2) 

.4[-  V)-(6+1)1i^  ‘^^6+1)] 

+6[-^h^“^\o)--|-  h^"^^{0)-6h^"^^(0)-h^"^^(0)+h^"'^^(6)] 

-4[-i^liLh^"^)(0)-ii^il-h^’^)(0)-(6-l)h^"^\o)-h^"^^(0)+h^'^)(6-l)] 


. h^" ^ 0 ) -  0) -( 6 -2 )h^  " 0 ) -h^  ~^\ 0 )+h^ 6 - 2) 

mg=h^  6+2)  -4h^  6+l)+6h^  ”‘^^(6)  -4h^  "‘^^(  6  -1  )+h^  6  -2) 


i 
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1 


For  the  other  three  submatrices  related  to  the  double-sided 


r  epr  e  s  entation. 


1  /  S  .  -8.  -js,^,  .  1  /  -St,  s.  ks  ds 

I'i-r  =  J  — (®  -2+e  )  e  ''  H(s)  -2+e  )“  - 

^  r  s  s 


H(a)f  2s  ,  s  ,  /  .  -8  ,  -2s\ 

^  e  -4e  +  6-4e  +  e  ■  e 


2s\  -(j+k)s  ds 

2TTJ 


=  0  all  j  >  1,  k  >  1 

The  elements  corresponding  to  j  =  1,  k  =  0,  and  j  =  0,  k  =  1  may¬ 
be  considered  as  belonging  to  one  of  the  other  submatrices 

t<  H(s)  f  2s  A  s,,  A  -s  .  -2s\  (j-fk)s  ds 
m~=  — V-M  e  -  4e  -{-6  -  4e  +  e  t  e  r  . 

Jk  p  ^4  I  J  2itj 

By  defining  6  =  j  -  k=  j  +  k  for  this  submatrix  and  noting  6  >  4, 
mg  can  be  computed  from  a  kno-wledge  of  ^kat  is, 

mg  =  h^-^)(6+2)-4h^"^\6+l)+6h^'^\6)-4h^''^^(6-l)+h^‘'^\6-2). 

This  completes  the  determination  of  the  M  matrix.  Due  to  the 
correlation  existing  between  |  and  triangular  matrices 

do  not  result  from  mj^  and  m^j^  as  in  Section  1,  4.  Apart  from  this, 
the  general  form  of  the  complete  matrix  is  very  similar  to  that 
derived  in  Part  1.4B. 

B.  The  H  Matrix 

To  obtain  the  H  matrix  it  is  necessary  to  pre-multiply 
the  M  matrix  by  G  where  G  is  the  relevant  Gram  Matrix.  Re¬ 
placement  of  H(s)  by  unity  in  the  previous  section  leads  to  a  tri¬ 
diagonal  matrix  with  elements: 

gg  =  2/3  6  =  0 

gg  =  l/6  6  =  ±  1 

gg=0  6<-l,  6>1 


1 


I 


Inversion  of  this  matrix  is  carried  out  in  Appendix  4,  from  which 


n 

-1  _  6  y 

’ij  “  5+T 


iktr  .  kdif 


k=l  2  + cos 


n+T 


n  =  2N  +  1 


-1  th 

where  is  the  element  in  the  ij  position  of  the  Inverted  Gram 

Matrix. 


I 

CO  « 


i  ; 
r 

is 


,) 


1 


1  2% 


1 


jH 


CO 
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Z.  CONTINUOUS  REPRESENTATIONS  OF  SYSTEM  MATRICES 
2.1  INTRODUCTION 

The  first  part  of  this  report  has  been  concerned  with 
the  reduction  of  a  known  system  transfer  function  to  matrix  form. 
This  part  considers  the  converse  problem  in  which  only  the  system 
matrix  is  known.  Since  this  matrix  contains  information  which 
relates  only  to  the  finite -dimensional  space  spanned  by  the  signal 
basis,  it  will  not  generally  be  possible  to  recover  the  original 
system  transfer  function  of  which  the  given  matrix  is  a  discrete 
representation.  In  many  cases,  however,  it  is  possible  to  find  a 
transfer  function  which  is  equivalent  to  that  of  the  actual  system 
in  the  sense  that  it  has  the  same  matrix  representation  with  respect 
to  the  relevant  signal  basis. 

Consider  the  case  in  which  this  equivalent  system 
transfer  function  is  described  by  the  infinite  sum  of  suitably  defined 
component  functions.  Formulas  derived  in  Part  1  of  this  report 
can  then  be  used  to  obtain  a  set  of  equations  relating  each  known 
matrix  element  to  the  unknown  coefficients  of  the  summation.  In 
general,  each  of  the  infinite  set  of  coefficients  will  contribute 
to  every  element  of  the  matrix,  and  the  relationship  may  not  be  a 
linear  one.  As  will  be  seen,  however,  there  is  a  considerable 
reduction  in  complexity  when  the  assumed  transfer  function  is 
represented  in  terms  of  an  infinite  set  of  the  same  component 
functions  as  describe  the  signal.  From  Equation  (1.  Z)  of  Part  1, 
we  have 

H  =  G"^  M 

where  H  is  now  the  known  system  matrix,  G  the  Gram  Matrix  for 
the  signal  basis  used,  and  M  consists  of  elements 

M  =  G  H 

so  M  is  readily  determined. 


't 


I 
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But 

where  k  =  1,  2,  ...  .N  is  the  signal  basis 

If  CO 

Hts)  =  ^  A.  4/.(s)  , 

i=l 

then 


2 

This  is  a  statement  of  the  N  equations  referred  to  above.  The 

evaluation  of  all  coefficients  A.  which  contribute  to  m,,  for  all  i 

1  ^  IK 

and  k  less  than  or  equal  to  N  determines  H(s),  giving  an  equivalent 
transfer  function  of  which  H  is  the  matrix  representation. 


CO 


12"  4 
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Z\} 


00 


2.2  EXPONENTIAL  EXPANSIONS 

Consider  a  signal  class  described  in  terms  of  forward 
exponential  functions 

""  7^  ’  ®k  ^  °  os  0,  k  =  1,  2, - N 

k 

when  the  s,  are- distinct, 
k 

Taking  the  diagonal  representation 
Hg  =  G  G‘* 

where  is  the  system  matrix  corresponding  to  the  basis  defined 
above,  we  have  N  equations  of  the  form 


h..  =  H(-s  .) 
JJ  J 


j  =  1,2,  ....N 


(2.1) 


If 


then 


00 

H(s)=  E 


•  j  S-8. 
1=1  1 


00 

E 


h..  = 

i  =  1  -s  '.  -s. 


A. 

1 

TT 


j  =  1,  2,  ...  .N 


J 


1=1  ^ 


(2.2) 


oo 

=  S  A  g 
i=  1  ^ 

If  H)  is  a  column  matrix  with  h..  in  the  jth  row  and  A)  is  a  column 

J  J 

matrix  with  elements  A^,  i  =  1,  2,  ...  .N,  then 


A)  =  G  H> 


(2.3) 


The  form  of  this  equation  suggests  that  for  the  assumed  ex¬ 
pansion,  the  coefficients  A^  are  those  which  provide  a  minimum 
least-squares  estimate  of  H(s).  Indeed  they  have  been  chosen  in 


such  a  way  that 


I 


-22- 

% 


A  >'fi  )!< 

H(-s  .)  =  H(-s  .) 

which  is  the  criterion  for  least-squares  coefficients  for  an  expan¬ 
sion  in  terms  of  a  preassigned  set  of  exponents,  s^,  j  =  1,  2,  ...  .N, 

An  alternative  approach  is  to  assume  a  rational  fraction  ex- 

A 

pansion  for  H(s),  that  is 


00 

S  b.'s^ 

Hf.)=  ‘ 


00 

Z  a.'s^ 
i=0  ^ 

and  to  make  use  of  the  known  properties  of  transfer  functions  of 
physical  systems.  Then,  from  Equation  (2.1),  we  have:- 

N/2-1 

S 


i=  0 

- W/z 


i 

b!( -8  .) 

1  J 


JJ 


i:  a!(-s.)^ 

i=0  ^  J 


n/2-1 


z  b.h..  -  a.  f-  ( -s  .f  =  (-S  .) 

^  1  JJ  iJ  '  j'  '  j' 


*.N/2 


(2.4) 


j  =  1,  2, . N 


where 


"i"  ^N/2 

and  N  is  even. 


and  b . 
1 


b! 

1 

^N/2 


This  set  of  N  equations  may  be  solved  for  a^,  and  b^,  i  =  0, 
4,2 . N/2  -  1. 

An  example  is  given  in  Appendix  5. 


A'- 


CO 


i 


1  24 


1 
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2.3  LAGUERRE  FUNCTION  EXPANSIONS 
Let 

.1,  /  \  T  /  \  X  /  (s4-a*)^"^ 

4j.(s)  =  L.(s)  =  +  (-a -a  )  ^ 

'  "  (8-a)' 


W 


O'r 


■  X> 


i  =  1,2, - N 


then,  from  Equation  (1.1):- 


h  -h  -(  a  ^  r  (-s+q/ '^(8+a*)^''’-{s+a*)^"^  ds 


hg  =  -(  -a  -a  ) 


*3/2  °° 


I  aJ 

i=i  ‘ 


(s-g)^  ^  ^  ds 
{s+a*)^-i+2  TFJ 


=  -(-a-aV/^J  [ 


(  -a  -a  ) 


r 

( s  -a  ' )  ( 8  -a  )  ( s 

,3A  [ 

■^6  ''^6+1  -] 

(a+g’V  (-a*-a)^J 

^  { 

■^6  ■  '^6+1} 

^*,1/2 

6  -1 

E  h 

i=0 

6+1 

"XT 


1  ds 


(3.1) 


(3.2) 


If  a  =  -1/2 


6-1 


A^  =  T.  h. 

^  i=0  ^ 


(3.3) 


To  establish  that  A  ,  6  =  1,  2,  ....  N  are  least-squares  coefficients  it 
6 

is  sufficient  to  show  that  least- square  coefficients  satisfy  Equation 
(3.1).  Consider  the  orthonormal  expansion 


I 


I 
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c» 


1  1 


H(8)  =  E  A!  .L,(s) 
1  =  1 

<lJh) 


H  ~  ^  H  ■  ^6+1 1 


.6-1  ,  .  j6 

(-s-aY  (-8-aV 


_/  ^  r  r  (-8+a)“  1  TT/n>  'i® 

^  ^  I  J  L  —■:  r„  _*.6+i  J  2ij 


{s-a)®‘^  (s-a)^  1u/„»  ds 

^  ■ :  ,  ^6+1 


(s+a  )  (B+a  )‘ 


2itj 


/  *.i/2  r 

=  -(  -a -a  )  J 


*\l/2  r  (9-a)®"^(8+a*-s+a)  da 


(8+a  ) 


H(8) 


da 

Zirj 


=  .(.a-a’V'^^  J  L.(-s)*H(s)I^{8)  ^ 


=  -(-a-a  )  '  hg 

Thus  a  choice  of  orthogonal  Laguerre  Functiona  aa  haaia  also  leads 
to  least-squares  coefficients. 

An  example  is  given  in  Appendix  6. 


OO 


I 


i.M 


■C.O 
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2.4  EXPANSIONS  IN  TERMS  OF  ELEMENTARY  HOLD  FUNCTIONS 

X  ^  A  1  /i  -(i-l)B 

Let  H(s)=  E  A.  — (1-e  )e 

.  •  IS 

1=1 


then,  from  Equation  (1.1) 


00 

S 

i=l 

-J 

r 

00 

Ij 

i=l 

1  Jp 

-2  , 

h  =  0 
0 

1  , 

h5  =  0 

z¥r 


r  3s  -  s ,  _  -s\  (6-i)8  ds 

|e  -3e  +3-e  j  >  — 


in  the  right  half  plane 


6-1=1,  1*5  =  J 


in  the  left  half  plane  nor  at  the  origin 
e"  ds 

■J  2iJ  = 


'■6+1 

TT" 


6-1=0,  hj  =  Aj,  I  (e2“-3e  =  ) 


2s  o  _  s,  ds 
2Trj 


"^6 

TT 


Vi  -  ^ 

{  -^6+1 

A,  =:  2 

0 

£  -..I 

E  (-1)®-'+^  h 
i  =  0 

(4.1) 

This  choice  of  basis  does  not  lead  to  least  squares  coefficients. 
An  example  is  given  in  Appendix  7. 


I 


CO 
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2.5.  EXPANSIONS  IN  TERMS  OF  ELEMENTARY  TRIANGULAR 
FUNCTIONS 


From  Equation  (5. 1)  of  Part  1: 


r  ® 


I  /  \  1  /  S  -  .  -8.  -is 

4Ji(s)  =  -^  (e  -2  +  e  )  e 
s 


From  Equation  (1.1):- 

i-i  ^  s  J 


i  s  6  -f  3,  mg  =  0 


i  ^  6  -  3,  mg  =  0 


i  =  6  -t  2, 
i  =  6  -t  1, 

i  =  6 

i  =  6  -  1 


since  integrand  has  then  no  poles  in  the  right 
half  plane. 

since  integrand  has  then  no  poles  in  the  left 
half  plane  nor  at  the  origin. 


”“6  ^  '^6+2, 


e°  ds  1 


8 


F  2TTj  5!  6  +  2 


1  ,  2s  ,  s.  ds  26 


"^6  =  ^6-nJ 

r  ® 
r  " 


2iTj  FT  6-H 
28  /  2s , , Sx  ds  _  66  . 


m_  =  A. 
o  o 


r  1  ,  4s  ,  3s,  -c  2s  8.  ds_26  . 
-ij  -F<^  +^5e  -20e 


.  p  1  ,  5s  /  4s ,  . 3s  2s ,  .  c  s,  ds  1  . 
'^6"^6-2j  -F^®  +^5e  -20e  ■fl5e  )^^Ag_. 


i  =  6  -  2, 

=  h  [^6-t-2  +  ^^^6+1  +  ^^^6  +  +  ^6-2] 


^6-i2  "  ■  ^6-2 
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Since  =  0, 

6  s  -  2,  successive  coefficients  are  readily  obtained: - 

=  5;' 

“-1 

a; 

A2=  5! 

m  -  2,6A. 

0  1 

A3=  5! 

m^  -  26A2  -  66A|^ 

A^=  5! 

m^  -  26A3  -  -  26A^ 

An  example  is  given  in  Appendix  8. 


f^T 
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CONCLUSION 

Part  1  of  this  report  deals  with  the  matrix  representation  of 
a  known  linear  system.  This  representation  gives  the  response 
of  the  system  which,  in  terms  of  the  components  of  the  input  signal, 
has  minimum  error  in  the  least-squares  sense,  and  input-output 
calculations  are  reduced  to  simple  operations  with  matrices.  Basic 
relationships  are  established  for  component  functions  which  are  not 
necessarily  orthogonal  and  a  brief  discussion  of  the  system  matrix 
for  exponential  basis  with  distinct  or  repeated  exponents  is  included. 
The  remainder  of  Part  1  deals  with  the  computation  of  matrix 
elements  for  signal  classes  for  which  two  commonly  used  time 
representations  are  appropriate.  In  one  case,  the  component  func¬ 
tions  are  not  orthogonal  and  the  required  Inverse  Gram  Matrix  of 
arbitrary  dimensionality  is  established  in  an  appendix. 

One  further  advantage  of  the  matrix  representation  is  that  it 
retains  only  those  system  characteristics  which  e  ffectively  modify 
the  signal.  This  means  that  given  a  system  matrix  it  is  not 
generally  possible  to  infer  the  overall  properties  of  the  actual 
system.  It  is  possible,  however,  to  determine  a  transfer  function 
which  has  the  same  matrix  representation .  If  this  Equivalent 
Transfer  Function  is  obtained  in  a  suitable  form,  instrumentation 
of  the  discrete  system  maybe  achieved  by  standard  synthesis 
techniques.  Further,  when  an  exponential  basis  is  used,  an 
approximation  to  the  actual  system  transfer  function  is  obtained 
with  expansion  coefficients  chosen  in  the  least  squares  sense. 

Part  Z  of  this  report  deals  with  the  derivation  of  Equivalent 
Transfer  Functions  which  are  represented  by  an  infinite  set  of 
the  same  components  as  describe  the  signal.  This  restriction 
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leads  to  simple  expressions  relating  matrix  elements  and  expansion 
coefficients  for  all  signal  bases  considered  in  Part  1.  Examples 
are  given  in  appendices.  It  should  be  noticed  that  the  computation 
of  these  coefficients  is  accompanied  by  a  rapidly-increasing  loss 
of  accuracy.  Considered  in  relation  to  the  accuracy  with  which 
the  matrix  elements  are  known,  this  factor  sets  an  upper  limit 
on  the  dimensionality  of  the  discrete  system  for  which  this 
representation  is  useful. 
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APPENDIX  I,  ORTHONORMALIZATION  AND  THE  GRAM  MATRIX 


'®ik“  by  definition 

For  an  orthonormal  basis 

<?i  I "  ®ik 


Let 


N 

E 

i  =  1 

‘^ikl^: 

N 

^  1 

E 

(4',  c. 

J  J 

-  fi 

—  u 

Ik  "" 

N  N 


i=l  j=l 


■ji  '"ik  ®ji 


1  =  C  G  C 

G  =  £-‘c-‘  =  [CC] 


This  gives  the  relationship  between  the  Gram  Matrix  and  the 
matrix  of  the  transformation  which  provides  an  orthonormal  basis. 
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APPENDIX  IL  CHOICE  OF  INTEGRATION  PATH  IN  COMPUTING 
INNER  PRODUCTS 

00  oo 

<f’|g)  =  f  f(t)’'' g(t)  dt  =  f  f(t)  g(t)  dt  ,  for  f(t)  real. 

8=0 


F(s)  =  J  f(t)  e"®*'  dt  ,  Yj  cr  V2>  ® 
-00 

-St 

G{s)  =  J  g{t)  e  dt  .  0  <  112 


(?|g>  =  f  f(tl[r'"  G(p)eP'|^].-dt 


C-JOO 


8=0 


C-JOO 


c+joo 

=  I  G(p)F{8-p)-& 


C-JOO 


Y^<  Re  (s-p)  <  Y2 


For  path:- 


i1l<c<ii2 


also, 


Y^  <  Re(-p)  <  Y2>  since  s  =  0 
that  is, 

-Vz  <  c  <  -Y^ 


12  4 


! 
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APPENDIX  III.  TWO  USEFUL  INTEGRALS 


J  HW  ,6,  ^  ^  I 


where  h(t)  is  the  Inverse  Laplace  Transform  of  H(s) 

f*  H{e)  6s  ds  ^  X  v-(-l)/x^xx 

J  ~r~  ®  lij  (6-X.)h'  '(X.)  +  J  h'  '(X.)d\ 


0  ‘'0 


where 


=  -6h^"^^{0)  +  -  h^'^^O) 


h<-^)(X)=  J  h(X)  d\ 


^  =  r  h(X) 
J  _4  2ttj  J- 


5  „6 


-_h(-')(X)i^  +r  ht-^Wi^dX 

=  -  L  h^'^^O.)  +  h^'^\x).^^g~P?-  +  F  h^‘^\x)(6-X)dX 

=  -L  h^'^^O)  -^-i-  h^"^^(0)-6h^‘^>(0)-h^'^\o)+h^'^)(6) 


S  ^  h^4-i)^0)  +  h(-4)(6) 
i=0 
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APPENDIXIV.  INVERSION  OF  THE  GRAM  MATRIX  FOR  A 
TRIANGULAR  BASIS 

From  Section  1.4B 


C'i 


'CC 


k 


Z/3 

l/6 

0 

0 

0 

0 

l/6 

z/3 

l/6 

0 

0 

0 

0 

1/6 

z/3 

l/6 

0 

0 

0 

0 

1/6 

z/3 

1/6 

0 

0 

0 

0 

1/6 

z/3 

1/6 

0 

0 

0 

0 

1/6 

z/3 

Let  F  =  6G  -  6\  1,  then 


c 

1 

0 

0 

0 

1 

c 

1 

0 

0 

0 

1 

C 

1 

0 

F  = 

0 

0 

1 

C 

1 

0 

0 

0 

1 

C 

.... 

..... 

.... 

•  •  • 

.  .  .  » 

where  c  =  4  -  6X. 

To  obtain  eigenvalues  it  is  required  to  solve  the  equation 


det  F  =  F  =0 
—  n 


n  =  ZN  +  1  where  N  is  the  dimensionality  of  the  basis 


F  =  c  F  .  -  F  , 
n  n-1  n-Z 


F  ,  .  =  cF  -  F  . 
n+1  n  n-1 


! 


-6- 


This  may  be  written: 


n 


0 


In  order  to  raise  this  matrix  to  the  nth  power  the  Cauchy 
Integral  Formula  will  be  used.  Discussion  of  the  general  conditions 
under  which  this  formula  may  be  applied  to  matrices  is  beyond  the 
scope  of  this  memorandum,  but  verification  of  the  result  obtained, 
is,  of  course,  possible. 

f(z) 


f(A)  = 


zl  -A 


R 


dz 

Zwj 


where  R  encloses  each  zero  of|zl  -  Ajand  no  singularities  of  f(z), 

-1 

n  dz 
^  Zirj 


n 

n 

_  r 

0 

1 _ 

1 

[l  0 

“  J 

R 

.-1 

z 

-  r 

n 

z 

z 

-1 

d 

J 

R 

z(z-c)+l 

1 

z-c 

2 

-1 

i 

c  =  a  +  a 


Let  a,  a  be  roots  of  z(z-c)  +  1=0 
-1 


- 

1 

n 

1 

r 

a 

c 

-1 

;l 

P 

a 

-1 

-n 

1 

_1 

o_ 

1 

P 

1 

P 

1 

a  -c_ 

P 

1 

1 

P 

_1 

P 

1 

1 

o 

l _ 

^n+l 

c 

-i' 

F 

n 

4; 

F 

n 

1 

r 

0 

r  T 

I 


JO 


a -a 


n+1  -n-1 

-n 

-  a 

-a  +  a 

n  -n 

n-1 

-n+1 

“  a 

-a  +a 

F 

r 

F 


since 


F  =1 
■^0 

For  F  =0 
n 


(a’^-a'”)(a+  a'^)  -  a“""  +  a'“'"  =  0 


n+1 

1 

"  n+1 
a 

-n-1 

a 

n 

a  -a 

n 

a 

-n 

a 

1  "" 

c  =  a 

+  a 

1 

n,  -n 
-a  +  a 


n-1  -n+1 

-a  +  a 


a  +a 
1 


-n 


-1.  n-1  .  -n-1 


that  is 


n+1  -n-1  n 

a  -  a.  =0 

2  n+  2  ■  j  2jf  ir 

a  =  1  = 

^  _  ^j21ir/2n+2  _  ^jlir/n+l 


1  =  1,  2,  ...  .n 


=  4  -  6X, 


^  gjiir/n+l  ^  g-jiir/n+l 


=  2  cos 


ilT 

n+T 


T  I  IT  1 

^  H+T  J 


\  =  T  [ 

or,  if  k  =  n  +  1-i 
X 

}f 

Coupled  Oscillator  theory  suggests  eigenvectors  of  the  form 

.  mkir 
u  =  sin  — r-j- 
m  n+1 


k  =  T  . N- 
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For  normalized  eigenvectors 


„  ,  .  mkir  . 

^  ^  II+T  "  ^ 

m=l 

It  can  be  shown^  that  this  requirement  is  satisfied  if 


that  is, 


u 


mk 


Vn+T 


sin 


mkir 

II+T 


■where 


A  = 


-1  rr 

-1  ,  - 

1 

=  U 

A 

U 

■^1 

0 

0 

. . . 

0 

^2 

0 

0 

0 

^3 

. . . 

. . . 

, . . 

.  .  . 

. . . 

. , , 

...  . 

. . . 

-1  _ 

N 

n 

Lj, 

S 

k=l 

E, 

i  =  l 

'^ik 

n 

=  s 

k=l 

^ik 

\k 

n+1 


n 


k=l 


,  ikir  .  kiir 
am  .  ■  sin— ^ 
n+1  n+1 


2  +  cos 


kir 

H+T 


This  gives  the  required  Inverse  Gram  Matrix. 


CO 


V 
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APPENDIXV.  EXPONENTIAL  EXPANSIONS.  EXAMPLE 


Hi  is  given  with  respect  to  an  orthogonalized  harmonic 
— (j) 

exponential  basis  by: - 


<* 


CO 


1,  166667 
-0.  942809 
0. 144338 
-0, 066667 


0. 

0.833333 
-0.  898146 
0.  282843 


0 

0 

0. 650000 
-0, 808290 


0 

0 

0 

0.533333 


(i)  From  Equation  (2.  3):- 


- 

- 

-1 

•  “ 

1/2 

1/3 

1/4 

1/5 

1 . 166667 

1/3 

1/4 

1/5 

1/6 

0.833333 

1/4 

1/5 

1/6 

1/7 

0. 650000 

1/5 

1/6 

1/7 

1/8  _ 

0. 533333 

200  -1200  2100  -1120~ 

' 1.166667  * 

-1200  8100  -15120  8400 

0.833333 

2100  -15120  29400  -16800 

0.650000 

-1120  8400  -16800  9800 

0.533333 

■  1. 0009  ■ 

0.9945 

0.0008 

-0.0064 

1.0009  .  1.9945^  0.0008  0.0064 

H(s)  ~  -g+l-  +  s+z"  +  -i+l - i+¥~ 
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(ii)  From  Equation  (2.4);- 

0.857445  b  -a  +  0.  857145  b.  -  a.  =  1 
o  o  11 

1.  200000  b  -  a  +  2.  400000  b.  -  2a.  =  4 
o  o  11 

1.538461  b  -  a  +  4.615383  b.  -  3a.  =  9 
o  o  11 

1. 875000  b  -  a.  +  7.  500000  b.  -  4a.  =  16 
o  1  11 

Solving,  we  obtain 

,  3.00008+3.9959  _  0.9969  ,  2.0032 

8  +2.9985+1.9980  s+u.vv^o  s+i.-^VVU 


CO 


1.2  4 
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APPENDIX  VI.  LAGUERRE  FUNCTION  EXPANSIONS  EXAMPLE 

is  given  with  respect  to  a  basis  consisting  of  orthogonal 
Laguerre  Functions  with  =  -  l/Z,  k  =  1,  2,  3,  4,  by 


1. 46667  0 

-0.76444  1.46667 

-0.34015  -0.76444 

-0.16458  -0.34015 


1. 46667 
-0. 76444 


1.46667 


From  Equation  (3.3) 

A^  =  1.46667 

A^  =  1. 46667  -  0.  76444  =  0. 70223 
A^  =  0. 70223  -  0. 34015  =  0. 36208 
A^  =  0.  36208  -  0. 16458  =  0. 1975 


H{8)~  1. 46667  L^(s)  +  0.  70223  ^2(8)  +  0.36028  1.3(8)  +  0. 1975  Lj^{s) 


where 


^  (8  +  1/2)^ 


i=  1,2,  3,4 


1 
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APPENDIX  VII.  EXPANSION  IN  TERMS  OF  ELEMENTARY  HOLD 
FUNCTIONS  EXAMPLE 


given  with  respect  to  the  "B"  basis  by:- 


0.93555  0 

0,77340  0.93555 

0.19759  0.77340 

0.06092  0.19759 


0  0 

0  0 

0.93555  0 
0.77340  0,93555 


then,  from  Equation  (4.1) 


A,  =  2h  =  1,87109 
1  o 

A^  =  Zh^  -  A^  =  0.32429 
A3  =  2h2  -  A^  =  0.  71947 

A^  =  2h3  -  A3  =  -0.59762 
A 

H(s)~  1,87109  B^(s)  -0.  32429B2(8)  +  0.  71947  83(8)  -  0.  59762  B4(s). 


* 


cc 


/ 


— •» 


I 


APPENDIX  VIII.  EXPANSIONS  IN  TERMS  QI'  EX..EMENTARY 
TRIANGULAR  FUNCTIONS.  EXAMPLE 


A  system  matrix  is  given  with  respect  to  a  basis  consisting  of 
elementary  triangular  functions  by 

To. 053333  -0.0766671 

^  Ml  .136667  n  I 


0. 356667 


Obtain  H(s)  in  the  form  S  A.  f  .(s)  and  show  that  the  matrix 
A  1=1  1  ^ 

corresponding  to  H(s)  wim  respect  to  the  given  basis  is 

’2/3  1/61  To. 053333  -0.076667' 

Mi=GHi  = 

[1/6  2/3 J  [1.136667  0.356667 


225000 


0.008333 


0.766667  0.225000 


From  Equation  (5.1) 


=  120  (0.008333)  =  0.999960 
A^  =  120  (0.  225000)  -  26(0.999960)  =  1.001040 
A3  =  120  (0.766667)  -  26(1.001040)  -  66(0.999960) 


=  0. 024360 


H(8)  ~  0.9996  \|/^(9)  +  1.  001040  1)2(8)  -  0.024360  t)  3(3) 


From  Section  1.5  of  Part  1:- 


i 
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=  0 

h^'^)(l)  =  1/5!  (0.996) 
h^'^)(2)  =  1/51  (30.99984) 
h<'^)(3)  =  1/5!  (210.0018) 

=  h^''*^(l)  =  0.008333 

m  =  h^‘'‘^(2)  -  4h^"^\l)  =  0.225000 
o 

m,  =  h^''‘^(3)  -  4h^"'^>{2)  +  6h<‘^)(l)  =  0.765018 
1 


2/3 

1/6' 

-1 

0,  225000 

0. 

o 

O 

OD 

UJ 

U) 
- » 

1/6 

2/3 

0.  765018 

0. 

225000 J 

1.6 

-0.4' 

‘O.  225000 

0. 

o 

o 

00 

-0.4 

1.6_ 

0.  765018 

0. 

225000 J 

0.05393  0.76667 

1,13403  0.35667 


which  may  be  compared  with 
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